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A quantum Monte Carlo simulation method has been developed and applied to study the critical behavior of a single Kondo 
impurity in a Luttinger liquid. This numerically exact method has no finite-size limitations and allows to simulate the whole 
temperature range. Focusing on the impurity magnetic susceptibility, we determine the scaling functions, in particular for 
temperatures well below the Kondo temperature. In the absence of elastic potential scattering, we find Fermi-liquid behavior 
for strong electron-electron interactions, Qc < 1/2, and anomalous power laws for 1/2 < Qc < 1, where is the correlation 
parameter of the Luttinger liquid. These findings resolve a recent controversy. If elastic potential scattering is present, we find 
a logarithmically divergent impurity susceptibility at gc < 1/2 which can be rationalized in terms of the two-channel Kondo 
model. 

PACS numbers: 71.10.Pm, 72.10.Fk, 72.15.Qm, 75.20.Hr 



I. INTRODUCTION 

Since its discovery, the Kondo problem- is one of the 
central topics in condensed-matter physicsE'Hlt describes 
a magnetic spin-i impurity embedded into a metal and 
may be the simplest example for the growth of an ef- 
fective coupling at low energies, resulting in a nonper- 
turbative ground state. For normal metals this ground 
state is found to be of Fermi- liquid type, where the cuiasi- 
particle wave functions simply acquire a phase shift. 13 The 
situation might change in one-dimensional (ID) systems 
which are known to exhibit non- Fermi liquid behavior for 
arbitrary Coulomb interactions. The fundamental the- 
ory of interacting ID metalS|-jii, the low-energy regime 
is the Luttinger liquid model. eTeI It is therefore of inter- 
est to understand the Kondo effect in a Luttinger liquid. 
An additional motivation arises from recent advances in 
nanofabrication whick now allow for controlled experi- 
ments on ID systemsJj In the future the question of how 
magnetic impurities behave when coupled to ID metals 
might be of |Grucial importance for experiments on quan- 
tum wiresJJ^B carbon nanotubes,cLDr for edge states in the 
fractional quantum Hall regime.EJ 

The Luttinger liquid model unifies the low-tem- 
perature physics of many microscopic lattice models for 
strongly correlated fermions, with only very few phe- 
nomenological parameters. In particular, one has the 
dimensionless Coulomb interaction strength parameters 
gc and gs for charge and spin sectors, respectively, and 
the charge- and spin-density velocities Vc and Vg. The 
crucial assumptions are the absence of lattice instabili- 
ties (like Umklapp scattering), the absence of electron- 
electron backscattering, and that the Coulomb interac- 
tion potential is screened by mobile charge carriers close 
to the ID metal. As a simple model for interacting 



fermions, the Luttinger liquid model is widely used to 
study the influence of electronic correlations on dynami- 
cal properties of ID metals, in particular in the presence 
of impurities. Tiie case of a spinless impurity is by now 
well understood.E2l If the impurity has internal degrees of 
freedom, the situation is more complicated and subject 
of this paper. 

A Kondo impurity coupled to a_Luttinger liquid was 
first considered by Lee und Toner .E3 Employing the per- 
turbative renormalization group they established how the 
Kondo temperature Tk depends on the exchange cou- 
pling constant J. This turns out to be a power-law de- 
pendence, while for normal metals Xf^ is an exponen- 
tial function of the coupling constant .□ The-same power 
law was found by Furusaki and Nagaosa.ll3 These au- 
thors derived the correct SU(2) invariant scaling equa- 
tions in the weak-coupling regime and tentatively ex- 
tended them to the strong-coupling regime, where a sta- 
ble strong-coupling fixed point was found for both an- 
tifcrromagnetic and ferromagnetic exchange couplings. 
This strong-coupling fixed point describes a many-body 
singlet formed by the impurity spin and the conduction 
electrons, similar to what happens in a normal metal. 
Moreover, Furusaki and Nagaosa made detailed predic- 
tions concerning the low-temperature critical properties 
of the impurity, e.g., the magnetic susceptibility, the heat 
capacity, and the conductance. These quantities were 
found to exhibit power-law behavior with interaction- 
dependent exponents. 

However, it remained unclear whether the extra- 
polation of the perturbative scaling equations into the 
strong-coupling regime is justified. Recent boundary 
conformal field theory (CFT) results by Frojdh and 
JohannessonEd allow only two possible scenarios. Either 
the system belongs to the Fermi liquid universality class 
or it indeed has the properties predicted by Furusaki and 
Nagaosa. CFT itself is however unable to unambiguously 
decide which universality class is ultimately realized for 
the Kondo problem in a Luttinger liquid. Several re- 
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cent papers seem to favor the local Fermi-liquid picture. 
Schiller and Ingersent have discussed a truncated bpt 
related model which exhibits Fermi-liquid behavior.ta 
In addition, according to the numerical density-matrix 
renormalization group (DMRG) calculation of Wang,t£l 
Fermi- liquid behavior holds for a spin-i impurity inter- 
acting with a ID Hubbard chain. Recently, Chen et 
al. deduced Fermi-liquid laws from the-^rity and spin- 
rotation symmetry of a related model.llZl Here we shall 
address and resolve this controversial issue. 

So far few studies have dealt with magnetic impuri- 
ties exhibiting elastic potential scattering in. addition to 
the conventional (Kondo) exchange coupling.HEj Why it 
should be considered at all becomes clear after the follow- 
ing discussion. If one starts out from the usual Ander- 
son model to describe a localized orbital interacting with 
conduction electrons, the natural generalization in ID 
would include JJoulomb interactions among the conduc- 
tion electrons .l3 For uncorrelated conduction electrons, 
one can then derive the usual Kondo exchange coupling 
in the local moment regime [which is realized for large 
on-site repulsion and a single-particle impurity level deep 
below the Fermi, energy] by applying the Schrieffer- Wolff 
transformation.B This transformation generates the ex- 
change coupling [see Eq. ( |2.6| ) below] and, in addition, 
an elastic potential scattering term. In the correlated 
case of interest here, this latter term may be crucial 
since elastic potential scattering is relevant in a Luttinger 
liquidO What can happen to a magnetic impurity in a 
Luttinger liquid in the presence of strong potentiaJ scat- 
tering was first studied by Fabrizio and Gogolin.E2l They 
predict that at the strong-coupling point for the Kondo 
effect discussed above elastic potential scattering is ir- 
relevant for rather weak repulsive Coulomb interaction, 
namely for 1/2 < < 1- In contrast, for strong enough 
interaction, gc < 1/2, the potential scattering breaks up 
the system into two independent chains. The magnetic 
impurity then interacts with two subsystemS|4f hannels) , 
and the two-channel Kondo picture emerges.EJ 

In this paper we present a path-integral quantum 
Monte Carlo (QMC) method allowing for the computa- 
tion of thermodynamic properties of a Kondo impurity 
in a Luttinger liquid. This method is numerically exact 
within the statistical error bars inherent to the MC tech- 
nique. The main advantages of our method are the ab- 
sence of any system-size restriction, contrary DMRG 
simulations or several QMC lattice algorithmsjH3 and the 
possibility of treating arbitrarily correlated conduction 
electrons. However, as is well known, QMC simulations 
of spin systems of tj en have to deal with the fundamen- 
tal sign problem.c3'cJ It is caused by sign alternations of 
the QMC weight function for different system configu- 
rations. This generally leads to a small signal-to-noise 
ratio and hence to numerical instabilities. The approach 
we shall discuss here is also plagued by a sign problem. 
However, our sampling technique moderates the prob- 
lem to an exjtent allowing us to treat sufficiently low 
temperatures .Ea Moreover, we have also applied filtering 



techniquescj which provide a general method to ease the 
mentioned difficulties. Yet the approach described below 
suffers only from a minor intrinsic sign problem, and the 
use of the filtering technique is not really crucial to ob- 
tain the results described below. Finally we mention that 
for a calculation of the Kondo screening cloud around 
the impurity, a simpler version of the present simulation 
method has been employed by Egger and Schoeller.EEl 

One might ask why we chose to develop a new algo- 
rithm for studying the Kondo effect even though the ex- 
ceptionally stable and widely ased QMC impurity algo- 
rithm due to Hirsch and FyeE3 is available. The rea- 
son is that we have to include the Coulomb interactions 
among the conduction electrons which are responsible for 
the Luttinger liquid state. In the Hirsch-Fye algorithm, 
one traces out the conduction electron degrees of free- 
dom away from the impurity and then updates only the 
arising fermion determinant. However, by construction, 
this procedure works only if the conduction electrons are 
in the Fermi liquid state. By employing the bosonization 
method, as is shown below, one can in fact follow a sim- 
ilar route as in Ref. ^ and trace out the now correlated 
conduction electrons away from the impurity. From a 
historical point of view, it may perhaps seem surprising 
that the QMC technique has not been a major tool in 
the resolution of the conventional Kondo problem where 
conduction electrons are described by a Fermi liquid. As 
realized in Ref. |2^, the main obstacle is the exponentially 
small Kondo temperature, which in turn requires the 
study of extremely low temperatures difhcult to achieve 
in path-integral MC calculations. For the Kondo effect 
in a Luttinger liquid, however, the Kondo temperature is 
much higher and QMC simulations become feasible even 
in the asymptotic low-temperature regime. 

The outline of this paper is as follows. In Sec. II we 
discuss the Luttinger liquid model with a Kondo impurity 
and describe our Monte Carlo algorithm in some detail. 
In Sec. Ill results for a Kondo impurity in the absence 
of elastic potential scattering are presented, and Sec. IV 
gives results in the presence of additional strong elastic 
potential scattering. Finally, some concluding remarks 
are offered in Sec. V. 



II. THEORY AND QUANTUM MONTE CARLO 
METHOD 

The low-energy properties of correlated ID systems are 
most convenieiijtly described in terms of the bosoniza- 
tion method.oia The spin-i electron field operator is ex- 
pressed in terms of spin and charge boson fields which 
obey the algebra (we put Ti = 1) 

[0,{x), ^j{x')] = ^'-S,,sgn{x - x') , (2.1) 

where i, j denote the charge (c) or spin (s) degrees of free- 
dom. The canonical momentum for the ipi phase field is 
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therefore Ili{x) = dx^Piix). The two kinds of phase fields 
are not independent but basically dual fields. Written in 
terms of the boson phase fields, the right- or left-moving 
[p = ±) component of the electron annihilation operator 
for spin a — ± takes the form 



2ttvf 



rjpa exp 



X exp ipkpx + ip\/Trj2[ipc{x) + (Tips{x)] . (2.2) 



The bandwidth cutoff is Uc, and we put Uc = vphp = I in 
what follows {vp is the Fermi velocity). A corresponding 
latti ce c onstant can then be defined a.s a — vp/uc- In 
Eq. (2^), we have also included real Majorana fcrmions 
rjpa- Their purpose is to ensure proper anticommutation 
relations between operators for different branches labeled 
by pa. Since only products rip„7j±p±a- will appear in the 
Hamiltonian, a convenient choice for these products is 
(see also Ref. ^ 



(2.3) 



where are the usual Pauli matrices. There is a simple 



way to realize why Eq. (2.2) holds. Keeping in mind that 
"Hpaflpa — 1 for all p and a, one can easily check that all 
products of the operators TipaVp' <t' 

defined in Eq. (|2j) 
fulfill the correct algebra required by a nticommutation 
relations for the rip^r. Actually, Eq. (2^) shows only one 
of several possibilities to choose representations of Ma- 
jorana fermion products. Of course, one can verify that 
the subsequent results do not depend on which one we 
choose. 

Under the conditions specified in the Introduction, the 
effective low-energy Hamiltonian for the clean electronic 
system takes the simple. Gaussian form of the bosonized 
Luttinger liquid model 
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dx [n^^ + gT\dxcp,f] 



(2.4) 



In a system with full (Galilean) translation invariance, 
the velocities Vc and Vs are required to fulfill Vi = vf/ Qi- 
We have assumed this relation in Eq. (2^), bearing, in 
mind that for lattice models it need not be fulfiUedfl A 
general rule of thumb for the dimensionless interaction 
strength parameter is 
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where U is the forward-scattering amplitude of the 
screened Coulomb interaction. In the important case 
of repulsive interactions, < 1. The spin parameter 
should be set to = 1 in order Ip respect the underly- 
ing spin isotropy of the electrons B This value is also the 
fixed point value of the renormalization group (RG) if 



one incorporates electron-electron backscattering. In the 
remainder, we shall put gg = \ and neglect backscatter- 
ing. Following the usual perturbative RG analysis, jthis 
could at most lead to weak logarithmic correctionscl to 
the power laws found below. 

Next we consider what happens once a single magnetic 
impurity is brought into the Luttinger liquid, say, at a; = 
0. We envision a spin-^ impurity characterized by the 

spin operator 5 = ^r, where r denotes the vector of Pauli 
matrice s [th is is not to be confused with the appearing 
in Eq. ( |2.3| )]. In terms of the conduction electron spin 
density operator 



s{x) 



E 



(2.5) 



and a point-like exchange coupling J, the standard con- 
tact contribution to the Hamiltonian reads with s± = 

Hi = Js{0)S = Js40)Sz + ^ (s+(0)5_ + s_(0)5+) . 

(2.6) 

We consider only antiferromagnetic value s J > in this 
paper. Usin g the bosonization formula (2.2), the spin 
density (2_^) of the conduction electrons reads 



Szix) 



1 



dxVs{x) 



Tz sin[2/ci?a; + v 27r ipc{x)] cos[v 27r ipsix)] 



s±(x) — — exp\±V2jT i6s(x)]\ ±iTy cos\\/2jt ips(x)] 
na I 

+ sin[2fci?x + V2T:ipcix)]j . (2.7) 

Here the Tj m atrices come from the Majorana fermion 
products (IJ). 

Now we can incorporate elastic potential scattering. 
For that purpose, we need the bozonized form of the total 
electron density operator. 



p{x) = ^y2/Trdxipcix) 



(2.8) 



H Tz cos[2kFX + v^TTipcix)] sinlV 2Tr (ps{x)] , 

na 

where we have omitted the background charge density 
2fci?/7r. The 2fci? component stems from terms mixing 
right- and left-moving particles, while the slow compo- 
nent ^ dx^c comes directly from the densities of right- 
and left-movers. There is also a Akp component in p(x) 
not specified in Eq. ( |2.8| ) which dominates in the regime 
gc < 1/3. Since in that limit the Coulomb interac- 
tions are extremely strong, any elastic potential scatter- 
ing will be highly relevant. Including a point-like scatter- 
ing potential of strength V, one obtains first a forward- 
scattering contribution if = Vy^2/7rdxipci0)- This can 
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simply be absorbed by a phase shift in the sm[2kFx] or 
cos[2fc_Fa;] factors and is therefore omitted in the sequel. 
We are then left with the important backscattering con- 
tribution 

Hv ^ T, cos\V2^ipJ0)] sin\V2^(ps(0)] ■ (2.9) 

nvp 

For numerical calculations, it is advantageous to em- 
ploy a unitarily transformed picture such that the Hamil- 
tonian becmjjies explicitly real- valued. This is achieved by 
choosinep3 



(2.10) 



C/ = exp[V27ri6'«(0)5^] , 

such that 

Us±SzrU^ — — \ ±iTy cos[V2tt If s{x)] 



7ra 

Tx sin[2fcFa; 



(2.11) 



The total transformed Hamiltonian H = UHU^ then 
reads 

1 



H = Ho + 



J 



/27r 



JSzdx(ps{0) 



(2.12) 



H (TxSxSm[V2Tripc{0)] + Tj^S'^ cos[\/27ri^s(0)] 

TTVp \ 



sin[v27r93c(0)] cos[v27r(^s(0)] 

2V 

TTVp 



Tz cos[v27r93c(0)] sin[v27r(ps(0)] 



where the transformation leads to a change in the forward 
scattering, 



J = J — 2tivp 



(2.13) 



The unitary transformation (2.10) also removes the Os{Q) 
field from the Hamiltonian (in fact, it is just constructed 
to remove this phas e factor). 

From Eq. ( 2.12| ) it is obvious that the Majorana 
fermions are dynamically constrained to follow the im- 
purity spin dynamics since 



Therefore we must have 

- ± 2S, . 



(2.14) 



(2.15) 



The only manifestation of the Majorana fermions is the 
overall sign, which we set equal to +1 in the following. 

One can simplify the total Hamiltonian using prop- 
erties of the products (g) Sk appearing in Eq. ( 2.12 ). 
Evaluated in the jS'^jT^) basis, we find from Eq. ( 2.15 ) 



.\S.Tz) = \5{S.,~S'z) ^ {S'ASx\S.) 



^-5{S,,^S'z)^^{S'^Sx\S, 



{S'zt'z\S.®tAS.Tz) = \5{S.,S',) = \{S',\S,) . (2.16) 



Therefore we can reduce the original Hamiltonian 
through the substitutions 



Sy 

Sz 



Sx 

1/2. 



(2.17) 



This leads from Eq. ( |2.12| ) to 
1 



H = H, 



+ 







J 



12-K 



JSzdx^siQ) 



2^^^ {25':, (cos[%/2^(Pc(0)] - cos[%/2^(p,(0)] 

cos[V2^Vc(0)] cos[V2^</3s(0)]} 
2V 



2S, sin[V27r¥3c(0)] sin[V27rv3s(0)] 

TTVp 



(2.18) 



For further convenience, we have shifted the field 
by -^/tt/S. This changes sin[V27r93c] to cos[\/27r(pc] and 
cos[\/27riy9c] to — sin[V27r(y9c]. 

We now proceed by integrating out all boson fields 
93j(x) for X 7^ for a given impurity spin path, as these 
represent just Gaussian integrations. The Euclidean ac- 
tion can then be expressed as an average over new fields 
{t denotes Euclidean time extending from i = to 
t = l3 = l/fcsT) 



(2.19) 



qj{t) = \J2'K^pj{x = 0, i) , 



with the constraint being enforced by Lagrange multiplier 
fields Aj(i). Since the spin and charge modes are only 
coupled through the terms ~ J and ~ in Eq. ( 2.18 ), 
the elimination of the degrees of freedoms can be car- 
ried out independently for j = c and j = s. As the 
computation for the charge part follows the same line of 
reasoning as for the spin part (it can be obtained by re- 
taining Qc factors and disregarding the ^ J term), we 
only discuss the elimination of the ips field in the follow- 
ing. Results for the c field are then recovered at the end. 



see Eq. (2.23) 



After a partial integration, we have to integrate out 
(Ps{x, t) from a problem characterized by the effective ac- 
tion 

^cff = ^ y dxdt [{dt^s)^ + vl^idxVsf] 



dxdt LPs {x, t)Sz {t)S' (x) 



+ i J dt Xs{t)[qsit) - V27rv3,(0,i)] . 

This can be achieved by solving the Euler-Lagrange equa- 
tion 



idf+v'pd'x)iPs(.X,t) 



-V2Tri 



J 



Xs{t)S{x)-t — S,{t)S'{x) 
Ztt 
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which is easily done in Fourier space, 

dk 



v. 



n—-~oo 



2tt 



iPsik,Un) , (2.20) 



with similar relations for other fields (aj„ = 27rn//3 are 
the Matsubara frequencies). Inserting the solution of the 
Euler-Lagrange equation for ips into ScS, one has 



+ (J/27r)2F;'(0,w„)5.K)5,(~a;„) . 

Here we have defined the boson propagatorsEl (j — c, s) 

exp{ikx) 



Fj{x,U!) = VF 



dk 



(2.21) 



exp(— |cja;/wj I) 



F'^{x,Lo) = {d''/dx^)F,{x,Lo) 



2'Kgj 



Six) 



exp(— |a;x/uj|) 



The (5(x)-contribution to F^'{x,u!) is irrelevant in our 
case, since it causes only a constant term ~ j dtSl{t) 
in the effective action. We therefore disregard it in the 
following. Finally, the Lagrange multiplier field can be 
integrated out by simple minimization, 



\s{ljJn) = ~i 



2.22) 



Collecting results, the effective action is found to read 



j — c,s n 

+ J dt cos[(?c(i)] cos[(7s(i)] 

/ dt2S,it)sm[q,{t)]sm[q,{t)] 

nvF J 



(2.23) 



+ 



Snf] 



with 5" formally given as J dtH'(t) with 
J 



H'it) 



2tivf 



-2S^ 



,(t)|cos[(7c(i)] - cos[(7,(i)]} . (2.24) 



After these preparations, we now proceed further and 
describe a quantum Monte Carlo (QMC) algorithm for 
this problem. Since the unitary transformation U given 
in Eq. ( ^.10 ) leads to the real- valued Hamiltonian (2.18), 



it is very convenient to employ this representation. We 
have focused on the impurity susceptibility 



X = 



/3 



dt{ss)Szm , 



(2.25) 



since knowledge of x at low temperatures is sufficient to 
answer the questions raised in the Introduction. Here 
the average is taken using Eq. ( 2.18| ). The impurity spin 
operator Sz does not change under the unitary trans- 
formation [/, so the expression ( 2.25| ) holds also in the 
transformed picture. 

The QMC simulation scheme starts out from the dis- 
creti zed im aginary-time path-integral repr esentati on fo r 
Eq. ( ^.25D using the effective action ( ^^2^ ) with ( ^.2^ . 
The imaginary-time slice is 5t — (3/N, where the Trotter 
number N should be large enough. In practice, one has to 
check empirically at the end that results converge upon 
increasing N. In the QMC simulations, a hard cutoff was 
chosen by keeping only Matsubara frequencies \ujn\ < t^c- 
The sampling of the qj fields is then most conveniently 
carried out directly using their Matsubara components. 
For the impurity spin variable, however, it is mandatory 
to use the time representation because one has a discrete 
variable Sj — 2Sz{tj) — ±1, where tj = jdt is the jth 
time slice. The action contribution S' is now determined 
as follows. From a Trotter breakup procedureEj valid at 
small enough St, we obtain the representation 



N 

exp(-y) = l[{Sj+i\eM-StH'{t,)]\Sj) , (2.26) 

where the spins obey periodic boundar y co nditions, 
Sn+i = Si. Using the matrix elements ( 2.16 ), we ob- 
tain (up to an irrelevant overall constant) 

AT 



exp{-S') = H {e/(*^) + 5,5,+ie-/fe)} , (2.27) 
where 

fit) = ^{cosbc(t)] - cosfc(t)]} . (2.28) 

The QMC sampling is then drawn from the weight 
function 



V ~ I exp(-^eff ) 



(2.29) 



where Scs is specified in Eq. ( p. 23 ) together with 

Eq. ( 2.27 ). Since exp(— 5") can be negative, the simu- 
lations have to face the sign problem.Ej For not exceed- 
ingly large J and low temperatures, however, the sign 
problem is not severe and the QMC algorithm described 
here can be applied to a wide region of the parameter 
space without instabilities. Denoting the sign of the MC 
weight as 



= sgn|exp(-S")| , 



(2.30) 



5 



the MC denominator will then be (^p). The severity p£ 
the sign problem is usually measured in terms of (^p).E^ 
One way to weaken thfi-.sign problem is to employ the 
Mak filtering techniques which can improve the stabil- 
ity of the algorithm by about 20 to 30%. For the results 
presented below, this technical trick was not necessary, 
and good statistics can be acquired even without a filter- 
ing method. 

Of particular interest is the value of the impurity 
susceptibility x which is given by the temperature- 
dependent expression 



6t {£.p SjSi 



(2.31) 



Here the Monte Carlo sampling over the configuration 
space spanned by the varia bles i qc(LUn), qs{oJn), Sj} is car- 
ried out using the weight (|2.29|). 



at least 5 full MC passes to ensure statistical indepen- 
dence. Our code performs at an average speed of 1 CPU 
hour per 5000 samples (separated by 5 MC passes) on 
an IBM RISC 6000/model 590 workstation at the lowest 
temperatures under consideration. Results reported here 
require typically 10® samples per data point. 



III. CRITICAL IMPURITY DYNAMICS 
WITHOUT POTENTIAL SCATTERING 

In this section we study the case without potential 
scattering (T^ = 0), with particular emphasis on the con- 
troversy about the low-temperature scaling. From our 
QMC data we observe that all x{T) curves for different 
coupling constants J but at a given Coulomb interaction 
strength gc can be mapped onto a single universal scaling 
curve. 
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FIG. 1. Low-temperature behavior of the impurity mag- 
netic susceptibility at (/c = 1/4 and various values of the cou- 
pling constant J. Notice the semi-logarithmic scales. Dashed 
curves represent guides to the eye only. Vertical bars give 
standard deviation error bars due to the MC sampling. 

For the results presented below, the average sign is 
(Cp) ~ but in practice stable simulations can be 
carried out even for (^p) « 0.01 at the expense of 
long central processing unit (CPU) times. The Monte 
Carlo trajectory was drawn from the standard Metropo- 
lis algorithm.cS We have used local updates of the phase 
fields at X = 0, i.e., of the Matsubara components qdoJu) 
and qs{ujn) for \LOn\ < lOc, and of the impurity spin tra- 
jectory Sz{tj) = Sj/2 — ±1/2. Typical discretization 
parameters [for J/2Trvp = 0.1 and = 0] required to en- 
sure convergence to the continuum limit of the discretized 
path integral are UcSt ~ 0.3. The acceptance ratios for 
local updates of the {Sj} variables are rather low for the 
parameter values considered below, typically of the or- 
der of 5%. Therefore data are accumulated only after 
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FIG. 2. Scaling curve for Fig. ^. Notice the 

semi-logarithmic scales. The dashed curve is a guide to the 
eye only. 

For instance, our raw data for t/c = 1/4 are shown in 
Fig. |l], and the scaling curve for x is depicted in Fig. ||. 
Apparently, there is a universal scaling function / such 
that 



TKxiT) = fiT/TK) 



(3.1) 



Using this matching procedure, the Kondo temperature 
can be determined straightforwardly. On the other hand, 
the value xo = x(^ = 0) is finite and can be used to 
define the Kondo temperature as well. Indeed, the zero- 
temperature magnetic susceptibility is of the order of the 
binding energy of the many-body singlet state formed by 
the impurity spin and the conduction electrons. There- 
fore we can fix Tk alternatively as 



(3.2) 
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which impHes /(O) = 1 from Eq. From the zero- 

temperature hmit of x(2^) [which can be obtained quite 
accurately by extrapolation of the data] we can then read 
off Tk. By means of either of these two prescriptions, as 
is shown in Fig. ^ for gc = 1 /4, one can indeed verify the 
dependence of the Kondo temperature on the coupling 
constant predicted in Rcfs. U2 and O, 



K 



D 



( J 



\2-nvp 



2/(l-9c) 



(3.3) 



where D is of the order of the bandwidth cutoff uJc- Si- 
multaneously, one gets the universal scaling curve for 
given interaction strength g^. 

Now we wish to address the low-temperature critical 
behavior. The low-temperature form (T <C Tk) of the 
impurity susceptibility ( |3.lD can exhibit only two POSn 
sibilities allowed from conformal field theory (CFT).til 
Either one has (i) Fermi-liquid behavior, 



/(t/Tk) = i-ci(r/rK)' + 



(3.4) 



or (ii) the anomalous exponents predicted by Furusaki 
and Nagaosa,E^ 



/(T/TA') = i-c2(r/rA') 



(3.5) 



where C\ and C2 are positive constants. Obviously, at 
gc = 1/2 one must see the behavior. This is a check for 
our numerics which is indeed passed nicely. The results 
(taking J/2'kvf = 0.1) are presented in Fig. ^. In the 
inset we have depicted the dependence of the deviation 
(Xo — x)Tk on the thermal scaling variable T/Tk- As 
one can see, the correct power law emerges, provided 
we are below the Kondo temperature. 




FIG. 3. Normalized Kondo temperature Tk as a function 
of the exchange coupling J for Qc — 1/4. The normalization 
has been chosen such that Tk ~ 1 for J/2ttvf = 0.1. The 
crosses represent the evaluation of Tk from Eq. (3.1) and the 
circles fro m E q . (3.2). The straight line has slopand the circles 
from Eq. (3.2). The str aigh t line has slope 8/3 and represents 
the prediction of Eq. (3.3). Notice the double-logarithmic 
scale. 
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FIG. 4. Scaling curve for the temperature dependence 
of the impurity susceptibility at gc = 1/2. The inset shows 
the same data for {xo — X)Tk as a function of T/Tk at low 
temperatures again. The straight line in the inset has slope 
2. Notice the double-logarithmic scales. The dotted curve is 
a guide to the eye only. 

The same critical behavior is found for gc — 1/4 as 
shown in Fig. ||. At low temperatures, the universal scal- 
ing curve displays a behavior. The data shown here 
were obtained for J jl'KVp = 0.1, but by virtue of scal- 
ing the same curve is found for other J as well. From 
these results one might be tempted to infer Fermi-liquid 
behavior for all interaction strength parameters gc- How- 
ever, we find the Furusaki-Nagaosa T^Ib- law as soon as 
gc > 1/2, as shown in Fig. ^ for gc = 3/4. The slope in 
the inset of Figjq is 4/3, in accordance with the exponent 
found in Ref. OF 
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T/Tk 

FIG. 5. Same as in Fig. | but for Qc = 1/4. The straight 
line in the inset has slope 2. 
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FIG. 7. Scaling curve as in Fig. ^ but for = 0.001. The 
straight line in the inset has slope 2. 



X 




T/Tk 

FIG. 6. Scaling curve as in Fig. | but for = 3/4. The 
straight line in the inset has slope 4/3. 

In addition, we have analyzed the case of extremely 
strong interactions. This formally corresponds to — > 0. 
Assuming that this limit is analytical, we can find the im- 
purity susceptibility behavior from a study oigc — 0.001, 
see Fig. |7[ Again, in accordance with our previous anal- 
ysis, we find a law for the low-temperature suscepti- 
bility. 



Let us now discuss these numerical results. Our sim- 
ulation data for the i mpur ity susceptibility at T <C Tk 
obey the scaling form (3.1) with 



f{T/TK) - 1 - c,iT/TKf - c^^TITkY's^ + • • • . (3.6) 

Hence there are two leading irrelevant scaling fields ,0 
one describing Fermi liquid (Ai ) and one describing non- 
Fermi liquid (A2) behavior, where the latter one corre- 
sponds to the Furusaki-Nagaosa prediction. For < 
1/2, the Fermi-liquid term is more important and leads 
to the observed behavior at low temperatures. In 
contrast, for 1/2 < < 1, the non- Fermi liquid behav- 
ior predicted in Ref. |l3| dominates. 

The scaling form ( ^.(|> i s consistent with the conformal- 
field theory analysis .ESE3 The operator O2 conjugate to 
the scaling field A2 is produced by a composite boundary 
operator in spin and charge sectors, while the operator 0\ 
comes from a composite operator given by the products 
of energy- momenta tensors in spin and charge sectors. As 
these are descendants of the identity operator, their con- 
tribution to the susceptibility becomes linear in Ai, i.e., 
ci '--^ Ai. In contrast, the Furusaki-Nagaosa term scaling 
like r^/S": is quadratic in the corresponding scaling field, 
C2 ~ A2. The amplitude of this contribution vanishes 
~ (1 — (7c) as gc — *■ 1, thereby reproducing the correct 
Fermi-liquid behavior of the conventional Kondo effect 
for uncorrelated electrons. Parenthetically, we note that 
the scaling field A2 also produces a subleading T^^^l9<' 
law in the impurity specific heat. 

To summarize, at g^ < 1/2, the Fermi liquid behav- 
ior will always dominate. However, at sufficiently low 
temperatures, the Furusaki-Nagaosa exponents can be 
observed at 1/2 < < 1- This finding is in conflict with 
the recent numerical DMRG study by Wang,E3 which re- 
ports Fermi-liquid behavior for a spin-i impurity coupled 
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to a Hubbard chain. Notice that the interaction param- 
eter gc for the ID Hubbard model away from half-fiUing 
is always within the bounds 1/2 < gc < IB Most likely 
the discrepancy is caused by finite-size effects due to the 
short chain len gths used in Ref. [l^. The more compli- 
cated outcome (3^) alsoj^ows that the simplified model 
by Schiller and Ingersentlla does not capture all essentials 
of the Kondo effect in a Luttinger liquid. 

Finally, let us discuss our data for extremely strong 
interactions, gc 0. For the clean case, it is well es- 
tablished that the Luttinger liquid model for gc ^ 
is equivalent to the low-energy sector of the ID Heisen- 
berg spin chain.cl Assuming that this reasoning carries 
over if a magnetic impurity is present, the scaling of 
the impurity susceptibility observed here (see Fig. for 
gc = 0.001) should also describe the susceptibility of a 
spin-i impurity interacting with a ID Heisenberg chain. 
One has to be careful to couple the impurity to just one 
site of the Heisenberg chain, otherwise an additional po- 
tential scattering contribution will be present (see Refs. 
|33|-p6| and Sec. IV). A different result was reported very 
recently by Liu,Ell namely a T^/^ scaling of the impurity 
susceptibility at low temperatures. Unfortunately, the 
reason for this discrepancy is not clear at the moment. 



IV. CRITICAL IMPURITY DYNAMICS WITH 
POTENTIAL SCATTERING 



X 




T/Tk 

FIG. 8. Impurity susceptibility at (jc = 1/4 in the presence 
of elastic potential scattering, 2V/tivf = 0.2, for exchange 
coupling J/2-RVF = 0.08. The data points are given as filled 
diamonds. For comparison, the U = data from Fig. ^ are 
shown as open circles. The solid line has slope I/ttF (see 
text), and the dotted curve gives has slope I/ ttF (see text), 
and the dotted curve gives the susceptibility ( |4.l[ ) of a free 
spin. The dashed curve is a guide to the eye only, and Tk is 
computed for V = 0. Notice the semi- logarithmic scales. 



In this section the influence of elastic potential scat- 
tering on the critical properties of a spin-i impurity in 
a Luttinger liquid will be discussed. As already men- 
tioned in the Introduction, for sufhciently strong interac- 
tion strength, gc < 1/2, and for strong enough potential 
scattering, the system is expected to displajj physics fa- 
miliar from the two-channel Kondo model.Ell 

In Fig. ^, data are shown for gc = 1/4. At high temper- 
atures, with or without elastic potential scattering, the 
Curie susceptibility of a free spin is always approached. 



Xfrcc(T) = /3/4 



(4.1) 



However, while the impurity susceptibility displays a 
crossover to the finite value xo — 1/Tk at zero temper- 
ature for vanishing elastic potential scattering strength 
(U = 0), the behavior is drastically different if potential 
scattering is present (here, ^V/ttvf — 0.2). The impurity 
susceptibility does not saturate but continues to increase 
without bound when lowering the temperature. Since we 
have semi-logarithmic scales in Fig. |^, our data are ac- 
curately fittedJpy the susceptibility of the two-channel 
Kondo modelB3 



X{T) 



ttF 



In 



(4.2) 



where F = Pjli^'^v],. This value for F (see Ref. |9[ 
indeed obtained from the slope of the solid line in Fig. ||, 



IS 



The corresponding results for gc ~ 3/4 are shown in 
Fig. ^. The logarithmically divergent behavior in the 
presence of potential scattering is not found anymore, 
and the low-temperature impurity susceptibility satu- 
rates at a finite value xX ■ Since Xo ^ XOi we expect from 
Eq. (3^) that all effects of potential scattering can be 
incorporated by a renormalization of the Kondo temper- 
ature Tk to smaller values. Within statistical error bars, 
the data for 2V/nvF = 0.3 shown in Fig. ^ can indeed be 
scaled onto the V = data, and the scaling function / 
holds even in the presence of elastic potential scattering. 
Clearly, this finding is in sharp contrast to the case of 
strong interactions, gc — 1/4, where potential scattering 
drastically changes the temperature dependence of the 
impurity susceptibility. 

From these data and the arguments of Ref. |l], we 
then expect two-channel Kondo behavior and hence a 
logarithmically divergent susceptibility for all gc < 1/2. 
An important special case of this general result is recov- 
ered for gc which corresponds to the ID Heisen- 
berg chain. Using numerical methods, bosonization a: 
CFT techniques, Eggert and AffleclicJ and Clarke et al 
have shown that a spin-i impurity in a Heisenberg chain 
exhibits a logarithmically divergent impurity susceptibil- 
ity. Due to the specific coupling of the impurity to the 
ID spin chain in these studies, an additional elastic po- 
tential scattering was present besides the usual Kondo 
exchange coupling term. 
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X 




T/Tk 

FIG. 9. Impurity susceptibility at = 3/4 in the presence 
of elastic potential scattering, 2V/ttvf = 0.3, for exchange 
coupling J/2nvF = 0.1. Data points are given as filled dia- 
monds. For comparison, the V = data from Fig. ^are shown 
as open circles. The dotted curve gives the susceptibility ( [l.l| ) 
of a free spin. The dashed curves are guides to the eye only, 
and Tk is computed for V — 0. Notice the semi-logarithmic 
scales. 



V. CONCLUSIONS 

In this paper the critical behavior of a spin-^ impurity 
in a correlated one-dimensional metal (Luttinger liquid) 
has been investigated numerically. To circumvent finite- 
size restrictions, we have developed and applied a quan- 
tum Monte Carlo algorithm which allows to determine 
any finite-temperature equilibrium quantity of interest. 
Here we have focused on the impurity susceptibility x 
with particular emphasis on the low-temperature behav- 
ior well below the Kondo temperature. Let us briefly 
summarize the main findings emerging from our numer- 
ically exact analysis. 

If elastic potential scattering is ignored, the impu- 
rity susceptibility shows the scaling behavior Tkx{T) — 
/(T/Tk) with a distinct universal scaling function / for 
each dimensionless interaction strength gc- It may be 
worth mentioning that scaling holds even outside the 
asymptotic low-temperature regime T <ti Tx- Within 
error bars, all data can be scaled onto universal scaling 
functions as long as T <C Wc, where Uc is the bandwidth. 
Matching x(T) curves for different J (but at a given g^) 
onto a scaling curve also yields the correct power-law de- 
pendence of the Kondo temperature, Tk ^ j2/(i-ff<:)^ 
which was first given in Ref. |l2|. 

We have then used our algorithm to determine the crit- 
ical behavior for T <C Tk ■ Generally one finds power laws 
x(T) ^ {T/Tk)"^ with some exponent ry. At = 1/4 and 



gc = 1/2, we find rj = 2, but at gc = 3/4, a different ex- 
ponent r] = 4/3 is obtained. Our data are consistent with 
the simultaneous existence of two leading irrelevant oper- 
ators, one describing Fermi-liquid behaiaor (77 = 2), the 
other describing the Furusaki-Nagaosall3 anomalous ex- 
ponent rj — 1/ gc- At gc < 1/2, the Fermi-liquid behavior 
is dominant, but at 1/2 < gc < 1, one can indeed observe 
the (7c-dependent ejmoiients. These findings resolve the 
recent controversytinia about the low-temperature criti- 
cality of the Kondo effect in a Luttinger liquid. 

We have also studied the effects of elastic potential 
scattering using our mpaerical approach. As predicted by 
Fabrizio and Gogolin,Eil for sufficiently strong Coulomb 
interaction strength, gc < 1/2, the impurity susceptibil- 
ity exhibits a logarithmically divergent behavior. The 
X ''^ ln(l/T) scaling is a manifestation of two-channel 
Kondo physics caused by the effectively open boundary 
at the impurity site. In contrast, for 1/2 < gc < 1, the 
susceptibility saturates to a finite value at zero tempera- 
ture and potential scattering does not modify the critical 
behavior. 

To conclude, we have numerically examined the criti- 
cal scaling properties of the Kondo effect in a Luttinger 
liquid. An interesting question which has not yet been 
studied in detail is related to universality in the pres- 
ence of potential scattering, e.g., the existence of uni- 
versal scaling functions for the impurity susceptibility. 
Future applications of our Monte Carlo algorithm might 
also deal with the case of more than one impurity, or with 
a systematic study of other quantities like the impurity 
specific heat. 
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